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Abstract 

Newtonian cosmological perturbation equations valid to full nonlinear order are well known in the literature. Assuming the ab¬ 
sence of the transverse-tracefree part of the metric, we present the general relativistic counterpart valid to full nonlinear order. The 
relativistic equations are presented without taking the slicing (temporal gauge) condition. The equations do have the proper New¬ 
tonian and first post-Newtonian limits. We also present the relativistic pressure correction terms in the Newtonian hydrodynamic 
equations. 


1. Introduction 

Cosmological perturbation theory is an important theoretical tool in interpreting cosmological observations like 
the two-dimensional temperature and polarization anisotropies of the cosmic microwave background radiation, the 
three-dimensional distribution and motions of galaxies, distorted images of galaxies due to gravitational lensing, 
etc. The cosmological perturbation equations are well known in the Newtonian context to fully nonlinear order IjJ], 
whereas the counterparts in Einstein’s gravity are known in linear lEi and low-order perturbation approximation 
101 • Here, we present a self-contained summary of the basic equations of recently formulated fully nonlinear and 
exact cosmological perturbation theory in Einstein’s gravity (Section[3]l. Comparisons are made with the Newtonian 
(Sections |2] and nil and the post-Newtonian equations (Section |5]l. We also present the Newtonian equations in the 
presence of relativistic pressure (Sectionj^. 

2. Newtonian cosmological perturbation theory 

Newtonian cosmological perturbation equations in the spatially homogeneous and isotropic background world 
model are l[ll] 


a__ i 

Q + 3-g = —V ■ (pv), 
a a 

a 1 1 1 _ 

V H- -V -H -V ■ Vv = -Vt/ - —Vp, 

a a a ag 

-tU - -4nG (g- g). 


( 1 ) 

( 2 ) 

( 3 ) 
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These are the mass conservation, the momentum conservation, and the Poisson’s equations, respectively; 'g, 'p, v, and 
U are the mass density, the pressure, the peculiar velocity, and the perturbed gravitational potential, respectively; a(t) 
is the cosmic scale factor. We decompose the mass density and pressure to the background and perturbed parts as 

^ - g + 6g, J’ - p + Sp. (4) 

Evolution of the background world model is described by equation (1211 1 properly derived in Einstein’s gravity. 

3. General relativistic cosmological perturbation theory 

We consider the scalar- and vector-type perturbations in aflat background with the metric convention fill 

ds^ = - (1 -I- 2a) c^dP' - 2xicdtdx' -i- (1 -t- 2(p) dijdx'dxP (5) 

where a, ip and^, are functions of spacetime with arbitrary amplitudes; index ofxi is raised and lowered by dg as the 
metric. We ignored the transverse-tracefree (TT) part of the metric which is interpreted as the gravitational waves to 
the linear perturbation order. The spatial part of the metric is simple because, in addition to ignoring the TT part, we 
already have taken the spatial gauge condition without losing any generality to the fully nonlinear order Hi. 

We consider a fluid without anisotropic stress. The energy momentum tensor is given as 


Tah^ gC^UaUh + p (gab + UaUb) , (6) 

where tildes indicate the covariant quantities; Mq is the normalized fluid four-vector; p includes the internal energy; in 
explicit presence of the internal energy we should replace 


p i H- -r-n , 


(7) 

where p in the right-hand-side is the rest-mass density ||3l. We introduce the following definitions of fluid three- 
velocities 


__ Vi 1 ^1 

M/ = a—, -Vi = Vi = — 

c y N 


{I+2ip) Vi - -Xi 
a 


(8) 


where y is the Lorentz factor 


r = 



v’^Vk 

1 

1 

c^(l + 2ip) 

/l 

/, l+2i^/r‘ Wvi Xt ) 


1 cHi+2<p) y 

N- \ c a{l+2(p)/\c a(\+2ifi)) 


(9) 


and N is related to the lapse function in (fT^ . The velocitiesTi and v/ are more physically motivated ones fst]: T/ is the 
fluid three-velocity measured by the Eulerian observer, and v/ is the coordinate three-velocity of fluid; the indices of 
Vi, 'vi and v/ are raised and lowered by dg. 

We can decompose andT/ into the scalar- and vector-type perturbations to the nonlinear order as ||3l 


I (v') — I — (^) 

Xi^cX,i+Xi\ Vi = -Va + v\', 


withy^*' *'. s 0 s v^' *'.. We assign dimensions as 


[fl] = [gah] = \^a\ = [a] = [ip] = Dk'] - [177c] = 1,^ 

[x']=L, \x]^T, [?/c]=L, [k\^T-\ [Tab]^[gc^]^[p], [G^^T-\ 


( 10 ) 


( 11 ) 


where k, the perturbed part of the trace of extrinsic curvature, will be introduced below. 

Here we present the complete set of fully nonlinear perturbation equations without taking the temporal gauge 13] . 
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The definition of k: 


a 


1 


The ADM energy constraint; 


3 (a^ 87rG__ Ac" 


2\a 


a 

+ -K + 


Nj N{l+2ip) 

c^Aip 1 T 


3^+3 + 


X <fi,k 


a a^(l + 2ip)^ 6 


1 +2(p 

4^G(e +J)(r"-1) + ^ 


3 


2aHl+2(p)^ 4 


( 12 ) 


(13) 


The ADM momentum constraint: 


2 c 

3 2a2jV(l+2<p) 


^Xi + ^x\ik I + (e + 3 I 


a^N{l+2ifi) 






(l+2^)2 

The trace of ADM propagation: 


-^^j+jWjj + YT^v, 


N 


(tV.; +Xi<4>'^ - -ii^ix ‘P,k 


:K,k, 


1 la 


cP' 


N \a 
1 


-377 - -^- 2 -^Me + 3 ^A + KP + -k + 2-K + 


Pan 


- + SttgIp + 




.1). 


N 

c 


a a^N(\ + 2ip) 


XK,i + c 


The tracefree ADM propagation: 


+ 3 — — K + 


cx 


\Ndt a a^N{\+2ip) )\a^N{\+2tp) 


a^N(\+2ip) 


pNj 

\ +2ip 


+ Pk)k\. 


\ +2ip 


X (tV,./ +XN'' - \^‘p^vA 


fl2(l + 2(p) 


1_ (v'V/ - yA ^ + 2 (v'V; - i<5 'a] 


1 + 2^ \ 


3 2 / 


TV 


= SttG 


a4TV2(l + 2(pp 
2 


Up’^Xlk-Xkjx'^A 


, p\\ y^j 1 ./_2 A 

r c2)[c2(l+2^t>) 3^1 2 

{)A''XkH^,j -x‘’’‘XN,k +Xkjx'‘<fi’‘ -XpkX'vA + Aipp -/xkV’‘v.j) 


1 + 2i;tJ 
„2 


fl2(l + 2ipp \ +2tp 

The covariant energy conservation: 

1 


{v'‘v,j - j + ^ + (fjN‘ - ^PyKi^ 


d 

— + 


dt a(l + 2tp) 

^ {NP)^k ^ NP‘P,k j_ 

a{\ + 2ip) a{\ + 2ipp y 
The covariant momentum conservation; 


(A/:^ + V)v.|s + (e + 7tJj„(32_,) 


_5 

dt a{\ + 2ip) 


r =0. 


1 

ay 


1 


_5 _ 

dt a(l + 2ip) 
1 


[nP + -/)' 


(ayv,) + v*V 


cXk 


a^(l + 2(p) I a 


/, M PNip,i 


'p-1 a(l + 2ip) 


g+ 4 (ar 


TV _ Pi 

gZ2. P’‘ ^2 


_5 

dt a(l + 2ip) 




? =0, 


( 14 ) 


( 15 ) 


( 16 ) 


( 17 ) 


( 18 ) 
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where k'j and N are the tracefree part of the extrinsic curvature and the lapse function, respectively, with 


" a^NHl+2(py 


1^/’-' {xij +Xj.i) - ^x' ,ix\j - Y7^ {xij +Xj,i) - ^x‘ 


(1 H - 2^)2 


x‘XiV’-’<fj + nX'x-'VjVj K N = Jl + 2a + 


x'^Xk 


a^(l + 2(f) 


(19) 


Apparently the basic set of perturbation equations in Einstein’s gravity in (fT^ - lfTSl) looks quite complicated com¬ 
pared with the Newtonian ones in ([T])-([3]l. In fact equations (fT2l i- (fT8l l are redundant; for example, the Einstein’s 
equations imply the conservation equations. More importantly, in the above equations, we have not taken the tempo¬ 
ral gauge condition yet. Equations (fTSl i- lfTST i are presented without taking the temporal gauge (hypersurface or slicing) 
condition. As the temporal gauge condition we can impose any one of the following conditions 


comoving gauge : 

<) 

III 

p 

zero-shear gauge ; 

o' 

III 

N 

uniform-curvature gauge ; 

III 

p 

uniform-expansion gauge ; 

k = Q, 

uniform-density gauge : 

III 

p 


or combinations of these to each perturbation order; we may call these the fundamental gauge conditions. With 
the imposition of any of these temporal gauge conditions the remaining perturbation variables are free from the 
remnant (spatial and temporal) gauge mode, and have unique gauge-invariant combinations. Thus, we can regard 
each perturbation variable in those gauges as the gauge-invariant one to nonlinear order iSIl. 

To the background order, (fT3l l. (fTSl l and (fTTT l. respectively, give 



where A is the cosmological constant. In the Newtonian limit we ignore p compared with gc^. 


4. Newtonian limit 

The infinite speed-of-light limit leads to the Newtonian equations |Q]. That is, as the Newtonian limit we consider 
the weak-gravity, the slow-motion, negligible pressure and internal energy compared with the energy density, and the 
small-scale (subhorizon) limits 


V Vk 


1 


O' « 1, « 1, p^QC, -r-Il^l, 




» 1 , 


( 22 ) 


where k the comoving wave-number with A = -k^\ H = ala-, in the presence of the cosmological constant A, we 
consider ~ 87rGp. We identify 


a- — \u, ip-2y^ V* = v, (23) 

where U and V correspond to the Newtonian and the post-Newtonian perturbed gravitational potentials, respectively 
fin ; equation (fTSl l gives ip - -a, thus V - U. We can show that equations ([T])-(l3l) follow from equations (fTTl i. (fTSl l. 
and (flSl l. respectively, in both the zero-shear gauge (^ s 0) and the uniform-expansion gauge (k = 0). 
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5. First post-Newtonian approximation 

To the first post-Newtonian (IPN) order our metric convention is fin 


ds^ - 




Pdt^ - -2aPicdtdx‘ + a^ll + -^2V I dijdx'dx-'. 


(24) 


where the index of P, is raised and lowered by dij. Similarly as in the metric of perturbation theory in equation (|5ll, in 
the spatial part of the metric we already have taken spatial gauge conditions without losing any generality, and have 
ignored the TT-type perturbation (13] ■ The dimensions are the following 


[U] = [V] = c^, [P,] = m = c\ 


(25) 


Comparing equation (O in the perturbation theory with equation (l24]l in the PN approach, to the IPN order we have 


Qr = —2 


y_ ‘ (y2_2®) 


1 1 / d -It, 

, k^--\ 3-U + 3V+-P\\, Xi^-aPi- 


(26) 


Using the identihcations between the two approaches made in equation (I26I) . we can derive the IPN equations 
from the perturbation equations in (fT^ - lfTSl) : v, is used as the fluid three-velocity in ||3]- Equation (fThl) to IPN order 
gives V - U. From equations (fTTl i. (fTsT l. (fTSi) and (fT4ll . respectively, we can show 


rjJ \ / >7 \ /,l 


p| -H -V ■ H- 21 ] H- n] H- (s- -I- -V ■ \\p 


dt a 


a a 


(27) 


1 


1 . 


- (avi) -I- -U / -H - — = — 

a a a a g 


-v^Uj + -(0-U^). + - iaPi)- + W {Pi,k - Pkj) 

n n ^ a 


1 


+ -(p 2 + 4 f; + n+^)^^-P,(^ + -v-V)(-y" + 3 t/)-v,:i(^ + -?-V)p 


P P,i 


g) g \dt a 


d 1 


-2 


_l/5 1. 


g\dt a 


(28) 


4t/ H- AnG(g-g) = -41 4 [2A<t) - 2USU + iaP' ,.)’l + 3U + 9-U + 6-U 

L \ J n n 


gv^ + ^ (?n - pu) -H ^ (p - p) |, 


-hSttG 
1 


(29) 

(30) 


where the right-hand-sides are IPN order, and we have recovered 11 explicitly; = v^v^. These are the same as 
cosmological IPN equations in fi0] . To IPN order, from equation (O we have 


(31) 


V, = (l-—t/lvi-H—P„ v; = (l-^|v, 


where Notice that in order to derive the IPN equations we have not imposed the temporal gauge condition; 

in the PN expansion Xi in equation (l26l l is already higher order. In the above equations we can still impose a gauge 
condition (Qj 


-P' : + nU + m-U = 0, (32) 

a ’ a 

where n and m can be arbitrary real numbers. The uniform-expansion gauge takes n - 2 - m, and the transverse-shear 
gauge takes n - 0 - m. We note that the comoving gauge, the uniform-curvature gauge and the uniform-density 
gauge are not available in the PN approximation. 
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6. Relativistic pressure corrections 

We can derive relativistic pressure correction terms in equations ([TJ-lO by relaxing the conditions p gc^ and 
n/c2 1 in equation (l22ll . We take the zero-shear gauge by setting x = 0. From equations (fTTl l. (fTsT l. and (flST l. 
respectively, we have 


P\ 

l~ p\ 

£>+ ^ + -V- 


\ j a 

A c2/ _ 


1 2 
—-V-Vp, 


a 


a 1 1 l/l_p, 

V H- -V H- - V ■ Vv- WU - - - ::^—r -V p + -^V , 

a a a g + p/c^ \a 


—jU- -AnG (q- g), 


(33) 

(34) 

(35) 


where p now includes the internal energy. We note that the pressure corrections in the above equations differ from 
the known ones in the literature The term in the right-hand-side of equation (l33T l was not known in the literature. 
Notice, in particular, the absence of pressure correction term in equation (ITST i. The presence of a pressure correction 
term, -127rG(p - p), in the right-hand-side of equation dlST l was often suggested in the literature |@,[^. 


7. Discussion 

The Newtonian equations in ©-([Si) are derived as the 1/c —» 0 limit of Einstein’s theory in the zero-shear gauge 
and the uniform-expansion gauge. The IPN equations in (l27]) - (l30l) are derived as a weak gravity limit keeping c^^- 
order terms compared with the Newtonian limit. Both the Newtonian and IPN equations are correctly recovered from 
our fully nonlinear and exact formulation of cosmological perturbation in Einstein’s gravity in equations (fT2l) - (fT8l) . In 
all (Newtonian, IPN, and Einstein’s) cases the equations are fully nonlinear and exact. In the nonlinear cosmological 
perturbation equations in the Einstein’s gravity and in the IPN approximation we can still impose one temporal gauge 
condition suggested in equations (l20l l and (l32l l. respectively. 

As we consider nonlinear perturbations, ignoring the TT part can be regarded as a serious assumption restricting 
the range of applications of our nonlinear perturbation formulation in Einstein’s gravity. Except for this shortcoming, 
as we have not decomposed the background parts, the basic set in equations (fT^ - (fT^ can also be regarded as an exact 
one; by contrast, except for the mass conservation equation, the equations in the Newtonian and IPN approximation 
are valid only to the perturbed variables in cosmology; the cosmological background equations are subtracted by using 
the equation derived in Einstein’s gravity 101. Meanwhile, by setting the background order quantities as a s 1 and 
g = 0 = p, equations ©-I© and equations ([2 ^-(l30l l are valid for the Newtonian and IPN hydrodynamic equations, 
respectively, in the Minkowski background 

Derivation of the equations in these proceedings with some applications will be presented in 11 [B. 
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